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In the constraint database model, spatial and spatio-temporal data are stored by boolean combi- 
nations of polynomial equalities and inequalities over the real numbers. The relational calculus 
augmented with polynomial constraints is the standard first-order query language for constraint 
databases. Although the expressive power of this query language has been studied extensively, 
the difficulty of the efficient evaluation of queries, usually involving some form of quantifier elimi- 
nation, has received considerably less attention. The inefficiency of existing quantifier-elimination 
software and the intrinsic difficulty of quantifier elimination have proven to be a bottle-neck for 
for real-world implementations of constraint database systems. 

In this paper, we focus on a particular query, called the alibi query, that was proposed in the 
context of moving object databases and that asks whether two moving objects whose positions are 
known at certain moments in time, could have possibly met, given certain speed constraints. This 
query can be seen as a constraint database query and its evaluation relies on the elimination of a 
block of three existential quantifiers. Implementations of general purpose elimination algorithms, 
such as provided by QEPCAD, Redlog and Mathematica, are in the specific case, for practical 
purposes, too slow in answering the alibi query and fail completely in the parametric case. 

The main contribution of this paper is an analytical solution to the parametric alibi query, 
which can be used to answer this query in the specific case in constant time. We also give an 
analytic solution to the alibi query at a fixed moment in time, which asks whether two moving 
objects that are known at discrete moments in time could have met at a given moment in time, 
given some speed constraints. 

The solutions we propose are based on geometric argumentation and they illustrate the fact 
that some practical problems require creative solutions, where at least in theory, existing systems 
could provide a solution. 

Categories and Subject Descriptors: H.2.3 [Database Management]: Languages — Query lan- 
guages; H.2.8 [Database Management]: Database Applications — Spatial databases and GIS; 
F.4.0 [Mathematical Logic and Formal Languages]: General 

General Terms: Theory, Mathematical Logic, Query Languages, Spatial databases 
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1. INTRODUCTION AND SUMMARY 

The framework of constraint databases was introduced in 1990 by Kanellakis, Ku- 
per and Revesz [Kanellakis et al. 1995] as an extension of the relational database 
model that allows the use of infinite, but first-order definable relations rather than 
just finite relations. It provides an elegant and powerful model for applications that 
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deal with infinite sets of points in some real affine space R", such as spatial and 
spatio-temporal databases [Paredaens et al. 1994]. In the setting of the constraint 
model, infinite relations in spatial or spatio-temporal databases are finitely repre- 
sented as boolean combinations of polynomial equalities and inequalities, which are 
interpreted over the real numbers. 

Various aspects of the constraint database model are well-studied by now. For 
an overview of research results we refer to [Paredaens et al. 2000] and the text- 
book [Revesz 2002] . The relational calculus augmented with polynomial constraints, 
or equivalently, first-order logic over the reals augmented with relation predicates to 
address the database relations . . . ,Rm, F0(-|-, x, <,0, . . . ,Rm) for short, 

is the standard first-order query language for constraint databases. The expressive 
power of first-order logic over the reals, as a constraint database query language, has 
been studied extensively [Paredaens et al. 2000]. It is well-known that first-order 
constraint queries can be effectively evaluated [Tarski 1951; Paredaens et al. 2000]. 
However, the difficulty of the efficient evaluation of first-order queries, usually in- 
volving some form of quantifier elimination, has been largely neglected [Heintz and 
Kuijpers 2004]. The existing constraint database systems or prototypes, such as 
Dedale and Disco [Paredaens et al. 2000, Chapters 17 and 18] are based on gen- 
eral purpose quantifier-elimination algorithms and are, in most cases, restricted to 
work with linear data, i.e., they use first-order logic over the reals without mul- 
tiplication [Paredaens et al. 2000, Part IV]. Of the general purpose elimination 
algorithms [Basu et al. 1996; Collins 1975; Grigor'ev and Vorobjov 1988; Heintz 
et al. 1990; Renegar 1992], some are now available in software packages such as 
QEPCAD [Hong 1990], Redlog [Sturm 2000] and Mathcmatica [Wolfram 2007]. 
But the intrinsic inefficiency of quantifier elimination and the inefficiency of its 
current implementations represent a bottle-neck for real-world implementations of 
constraint database systems [Heintz and Kuijpers 2004]. 

In this paper, we focus on a case study of quantifier elimination in constraint 
databases. Our example is the alibi query in moving object databases, which was 
introduced and studied in the area of geographic information systems (CIS) [Pfoser 
and Jensen 1999; Egenhofer 2003; Hornsby and Egenhofer 2002; Miller 2005]. This 
query can be expressed in the constraint database formalism and, at least in the- 
ory could be answered, both in the specific and the parametric case, by existing 
implementations of quantifier elimination over the reals. The evaluation of the al- 
ibi query adds up to the elimination of a block of three existential quantifiers. It 
turns out that packages such as QEPCAD [Hong 1990], Redlog [Sturm 2000] and 
Mathcmatica [Wolfram 2007], can only solve the alibi query in specific cases, with 
a running time that is not acceptable for moving object database users. In the 
parametric case, these quantifier-elimination implementations fail miserably. The 
main contribution of this paper is a theoretic and practical solution to the alibi 
query in the parametric case (and thus the specific cases). 

The research on spatial databases, which started in the 1980s from work in geo- 
graphic information systems, was extended in the second half of the 1990s to deal 
with spatio-temporal data. In this field, one particular line of research, started by 
Wolfson, concentrates on moving object databases (MODs) [Giiting and Schneider 
2005; Wolfson 2002], a field in which several data models and query languages have 
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been proposed to deal with moving objects whose position is recorded at discrete 
moments in time. Some of these models are geared towards handling uncertainty 
that may come from various sources (measurements of locations, interpolation. ...) 
and several query formalisms have been proposed [Su et al. 2001; Geerts 2004; Kui- 
jpers and Othman 2007]. For an overview of models and techniques for MODs, we 
refer to the book by Giiting and Schneider [Giiting and Schneider 2005] . 

In this paper, we focus on the trajectories that are produced by moving objects 
and which are stored in a database as a collection of tuples (ti,Xi,yj), i = 0, ...,N, 
i.e., as a finite sample of time-stamped locations in the plane. These samples may 
have been obtained by GPS-measurcments or from other location aware devices. 

One particular model for the management of the uncertainty of the moving ob- 
ject's position in between sample points is provided by the bead model. In this 
model, it is assumed that besides the time-stamped locations of the object also 
some background knowledge, in particular a (e.g., physically or law imposed) speed 
limitation at location (xj, y.i) is known. The bead between two consecutive sample 
points is defined as the collection of time-space points where the moving objects can 
have passed, given the speed limitation (see Figure 1 for an illustration). The chain 
of beads connecting consecutive trajectory sample points is called a lifeline neck- 
lace [Egenhofer 2003] . Whereas beads were already conceptually known in the time 
geography of Hagcrstrand in the 1970s [Hagerstrand 1970], they were introduced in 
the area of GIS by Pfoser [Pfoser and Jensen 1999] and later studied by Egenhofer 
and Hornsby [Egenhofer 2003; Hornsby and Egenhofer 2002], and Miller [Miller 
2005], and in a query language context by the present authors [Kuijpers and Oth- 
man 2007]. 

In this setting, a query of particular interest that has been studied, mainly by 
Egenhofer and Hornsby [Egenhofer 2003; Hornsby and Egenhofer 2002], is the alibi 
query. This boolean query asks whether two moving objects, that are given by 
samples of time-space points and speed limitations, can have physically met. This 
question adds up to deciding whether the necklaces of beads of these moving objects 
intersect or not. This problem can be considered solved in practice, when we can 
efficiently decide whether two beads intersect. 

Although approximate solutions to this problem have been proposed [Egenhofer 
2003], also an exact solution is possible. We show that the alibi query can be 
formulated in the constraint database model by means of a first-order query. Ex- 
periments with software packages such as QEPCAD [Hong 1990], Redlog [Sturm 
2000] and Mathematica [Wolfram 2007] on a variety of beads show that deciding 
if two concrete beads intersect can be computed on average in 2 minutes (running 
Windows XP Pro, SP2, with a Intel Pentium M, 1.73GHz, 1GB RAM). This means 
that evaluating the alibi query on the lifeline necklaces of two moving objects that 
each consist of 100 beads would take around 100 x 100 x 2 minutes, which is almost 
two weeks, or, if we work with ordered time intervals and first test on overlapping 
time intervals, (100 + 100) x 2 minutes, which is almost 7 hours. Clearly, such an 
amount of time is unacceptable. 

Another solution within the range of constraint databases is to find a formula, in 
which the apexes and limit speeds of two beads appear as parameters, that para- 
metrically expresses that two beads intersect. We call this problem the parametric 
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alibi query. A quantifier-free formula for this parametric version could, in theory, 
also be obtained by eliminating one block of three existential quantifiers from a 
formula with 17 variables using existing quantifier-elimination packages. We have 
attempted this approach using Mathematica, Redlog and QEPCAD, but after sev- 
eral days of running, with the configuration described above, we have interrupted 
the computation, without successful outcome. Clearly, the eliminating a block of 
three existential quantifiers from a formula in 17 variables is beyond the existing 
quantifier-elimination implementations. In fact, it is well known that these imple- 
mentations fail on complicated, higher-dimensional problems. The benefit of having 
a quantifier-free first-order formula that expresses that two beads intersect is that 
the alibi query on two beads can be answered in constant time. The problem of 
deciding whether two lifeline necklaces intersect can then be done in time propor- 
tional to the product (or the sum, if wc first test on overlapping time intervals) of 
the lengths of the two necklaces of beads. 

The main contribution of this paper is the description of an analytic solution 
to the alibi query. Wc give a quantifier-free formula, that contains square roots, 
however, and that expresses the (non-)cmptiness of the intersection of two para- 
metrically given beads. Although, in a strict sense, this formula cannot be seen 
as quantifier- free first-order formula (due to the roots), it still gives the above 
mentioned complexity benefits. Also, this formula with square roots can easily be 
turned into a quantifier-free formula of similar length. At the basis of our solution 
is a geometric theorem that describes three exclusive cases in which beads can in- 
tersect. These three cases can then be transformed into an analytic solution that 
can be used to answer the alibi query on the lifeline necklaces of two moving objects 
in less than a minute. This provides a practical solution to the alibi query. 

To back up our claim that the execution time of our method requires milliseconds 
or less we implemented this in Mathematica and compared it to using traditional 
quantifier elimination to decide this query. We have included this implementation 
in the Appendix and used it to perform numerous experiments which only confirm 
our claims. 

We give another example of a problem where common sense prevails over the 
existing implementations of general quantifier elimination methods. This problem 
is the alibi query at a fixed moment in time, which asks whether two moving objects 
that are known at discrete moments in time could have met at a given moment 
in time. This problem can be translated in deciding whether four disks in the 
two-dimsnional plane have a non-empty intersection. Again, this problem can be 
formulated in the context of the contraint model and adds up to the elimination 
of a block of two existential quantifiers. Also for this problem we provide an exact 
solution in terms of a quantifier-free formula. 

This paper is organized as follows. In Section 2, wc describe a model for trajectory 
(or moving object) databases with uncertainty using beads. In Section 3, wc discuss 
the alibi query. The geometry of beads is discussed in Section 4. An analytic 
solution to this query is given in Section 5 and experimental results in Section 6 of 
our implementation that can be found in the Appendix. The alibi query at a fixed 
moment in time is solved in Section 7. 
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2. A MODEL FOR MOVING OBJECT DATA WITH UNCERTAINTY 

In this paper, we consider moving objects in the two-dimensional (x, ?;)-space 
and describe their movement in the (i, x, y)-space R x R^, where t is time (we 
denote the set of the real numbers by R). 

In this section, we define trajectories, trajectory samples, beads and trajectory 
(sample) databases. Although it is more traditional to speak about moving object 
databases, wc use the term trajectory databases to emphasize that wc manage the 
trajectories produced by moving objects. 

2.1 Trajectories and trajectory samples 

Moving objects, which we assume to be points, produce a special kind of curves, 
which are parameterized by time and which wc call trajectories. 

Definition 1. A trajectory T is the graph of a mapping / C R ^ R2 : t 1-^ 
a{t) = {ax{t), a.y{t)), i.e., 

T = {(i, a^it), ay{t)) G R x R^ | t e /}, 

where / is the time domain of T . □ 

In practice, trajectories are only known at discrete moments in time. This partial 
knowledge of trajectories is formalized in the following definition. If we want to 
stress that some t, x, y- values (or other values) are constants, we will use sans serif 
characters. 

Definition 2. A trajectory sample is a finite set of time-space points {(to,Xo,yo), 
(ti, Xi, Yi), (tAT, xjv, yw)}, on which the order on time, tg < ti < • • • < tyv, induces 
a natural order. □ 

For practical purposes, we may assume that the (t^, x^, yi)-tuples of a trajectory 
sample contain rational values. 

A trajectory T, which contains a trajectory sample {(to,Xo,yo), (ti,Xi,yi), 
(tAr,XAr,yAr)}, i.e., (ti , (tj ) , (ti ) ) = (ti,Xi,y,) for i = 0, ...,iV, is called a geospa- 
tial lifeline for this trajectory sample [Egenhofcr 2003]. A common example of a 
lifeline, is the reconstruction of a trajectory from a trajectory samples by linear 
interpolation [Giiting and Schneider 2005]. 

2.2 Modeling uncertainty with beads 

Often, in practical applications, more is known about trajectories than merely some 
sample points (ti,Xi,yi). For instance, background knowledge like a physically or 
law imposed speed limitation Vi at location (xi, y^) might be available. Such a speed 
limit might even depend on ti. The speed limits that hold between two consecutive 
sample points can be used to model the uncertainty of a moving object's location 
between sample points. 

More specifically, we know that at a time t, U < t < U+i, the object's distance 
to (xi,yi) is at most Vi(t — t^) and its distance to (xi+i,yi+i) is at most Vi(ti+i ~t). 
The spatial location of the object is therefore somewhere in the intersection of the 
disc with center (xi,yi) and radius Vi(t — t^) and the disc with center (xi+i,yi+i) 
and radius Vi(ti+i — t). The geometric location of these points is referred to as 
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a bead [Pfoscr and Jensen 1999; Egenhofer 2003] and defined, for general points 
p = (tp, Xp, Up) and q = (tq, Xq, yq) and speed limit Umax as follows. 




Fig. 1. A bead and a lifeline necklace. 

Definition 3. The head with origin p = {tp,Xp,yp), destination q = {tq,Xq,yq), 
with tp < tq, and maximal speed Umax > is the set of all points {t, x,y) E Rx 
that satisfy the following constraint formula^ 

'^B(t,x,y,tp,Xp,yp,tq,Xq,yq,Vnii,^) 1= {x - Xpf + {y - ypf < {t-tpfvl^^^ 

A (.T ~ Xqf + {y- yqf < {tq - t)^ V^^^ A tp<t<tq. 
We denote this set by B(p, q, Umax) or B{tp, Xp, yp, tq, Xq,yq, Umax)- □ 

In the formula ■^s{t,x,y,tp,Xp,yp,tq,Xq,yq,Vm&^), we consider tp, Xp,yp,tq, Xq, 
yq, Umax to be parameters, whereas t, x, y are considered variables defining the subset 
of R X R2. 

Figure 2 illustrates the notion of bead in time-space. Whereas a continuous curve 
connecting the sample points of a trajectory sample was called a gcospatial lifeline, 
a chain of beads connecting succeeding trajectory sample points is called a lifeline 
necklace [Egenhofer 2003]. 

2.3 Trajectory databases 

We assume the existence of an infinite set Labels = {a, b, ai, bi, a2, b2, ...} of 
trajectory labels, that serve to identify individual trajectory samples. We now define 
the notion of trajectory database. 



^Later on, this type of formula's will be roferod to as F0(+, X , <, 0, l)-formulas. 
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Definition 4. A trajectory (sample) database is a finite set of tuples {ai,tij,Xi_j, 
yij-,Vij), with i ~ l,...,r and j — 0^...,Ni, such that a.; G Labels cannot appear 
twice in combination with the same t-value, such that {(t^.o, Xi^Oi Yi, 0)1 (ti,i, x^^i, yi_i), 
{ti^Ni,'X-i,Ni,yi,Ni)} is a trajectory sample for each i = 1, ...,r and such that the 
Vj,j > for each i = 1, ...,r and j = 0, Nj. □ 

3. TRAJECTORY QUERIES AND THE ALIBI QUERY 

In this section, we define the notion of trajectory database query, we show how 
constraint database languages can be used to query trajectories and we define the 
alibi query and the parametric alibi query. 

3.1 Trajectory queries 

A trajectory database query has been defined as a partial computable function from 
trajectory databases to trajectory databases [Kuijpers and Othman 2007]. Often, 
we are also interested in queries that express a property, i.e., in boolean queries. 
More formally, we can say that a boolean trajectory database query is a partial 
computable function from trajectory databases to {True, False}. 

When we say that a function is computable, this is with respect to some fixed 
encoding of the trajectory databases (e.g., rational numbers arc represented as pairs 
of natural numbers in bit representation). 

3.2 A constraint-based query language 

Several languages have been proposed to express queries on moving object data and 
trajectory databases (see [Giiting and Schneider 2005] and references therein). One 
particular language for querying trajectory data, that was recently studied in detail 
by the present authors, is provided by the formalism of constraint databases. This 
query language is a first-order logic which extends first-order logic over the real 
numbers with a predicate S to address the input trajectory database. We denote 
this logic by F0(-|~, x, <, 0, 1, S) and define it as follows. 

Definition 5. The language F0(+, x,<, 0,1,5') is a two-sorted logic with label 
variables a,b,c, ... (possibly with subscripts) that refer to trajectory labels and real 
variables x, y, z, v, ... (possibly with subscripts) that refer to real numbers. The 
atomic formulas of F0(-|-, x , <, 0, 1, 5) are 

— P{xi, ...,Xn) > 0, where P is a polynomial with integer coefficients in the real 
variables xi, a;„; 

— a = b; and 

— S{a, t, X, y, v) {S is a 5-ary predicate). 

The formulas of F0(+, x , <, 0, 1, S*) arc built from the atomic formulas using the 
logical connectives A, V, ^, ... and quantification over the two types of variables: 3a;, 
Vcc and 3a, Va. □ 

The label variables are assumed to range over the labels occurring in the input 
trajectory database and the real variables are assumed to range over R. The formula 
S{a,t,x,y,v) expresses that a tuple {a,t,x,y,v) belongs to the input trajectory 
database. The interpretation of the other formulas is standard. 
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For example, the F0(+, x, <, 0, 1, S')-sentcnce 

3a3b{^{a = fe) A \/t\fx\fy\fvS{a, t, x, y, v) <-> S{b, t, x, y, v)) 

expresses the boolean trajectory query that says that there are two identical tra- 
jectories in the input database with different labels. 

When we instantiate the free variables in a F0(-|-, x , <, 0, 1, S')-formula (p{a, b, 
t,x,y,...) by concrete values a, b, t, x, y, ... we write (/^[a, b, t, x, y, ...] for the 
formula we obtain. 

3.3 The alibi query 

The alibi query is the boolean query which asks whether two moving objects, say 
with labels a and a', that are available as samples in a trajectory database, can 
have physically met. Since the possible positions of these moving objects are, in 
between sample points, given by beads, the alibi query asks to decide if the two 
lifeline necklaces of a and a' intersect or not. 

More concretely, if the trajectory a is given in the trajectory database by the 
tuples (a,to,xo,yo, vo), (a,tAr,XAr,yAr, Vat) and the trajectory a' by the tuples 
(a',tQ,Xg,yg, Vg), (a',t^^,x'^j,y^j, v^j), then a has an alibi for not meeting a' if 
ioT all i, < i < N ~ 1 and aU j, 0<j<M~ 1, 

B(tj,x,,y,,tj+i,x,;+i,yj+i,v,) n B(t^ , x^-, y^-, t^-+i, x^+i, y^-+i, v^-) = 0. (f) 

We remark that the alibi query can be expressed by a formula in the logic 
F0(+, X , <, 0, 1, S), which we know give. To start, we denote the subformula 

5'(a,ti,xi,yi, vi) A S'(a, ^2, 2:2, 2/2, ^2) A 

yt3\/x3yy3yv3{S{a,t3,X3,y3,V3) -n{ti < A^s < ^2)), 

that expresses that (ti,a;i,?/i) and (t2,X2,y2) arc consecutive sample points on the 
trajectory a by cr(a, ti, xi, wi, i2, 2^2, 2/2, ■^2)- 

The alibi query on a and a' is then expressed as iy9aiibi[a, a'] = 

-^3ti3xi3yi3vi3t23x23y2^V23t[3x[3y[3v[3t'23x'23y2M 

{a{a,ti,xi,yi,vi,t2,X2,y2,V2) A a{a' ,t[,x[,y[,v[,t'2,X2,y2,v'2) A 
3t3x3y{ti < i < t2 A t'l < t < t'a A 
(x ~ xif + {y- yif <{t- hfvl A {x - X2f + {y - 1/2)' < (^2 - tfvl A 
(x - x',f + {y-y[f<{t- t\fvf A (x - xi)^ + (y - y',f < {t', - tfy'^)). 

It is well-known that F0(+, x , <, 0, 1, iS')-expressible queries can be evaluated 
effectively on arbitrary trajectory database inputs [Paredaens et al. 2000; Kuijpers 
and Othman 2007]. Briefly explained, this evaluation can be performed by (1) 
replacing the occurrences of S{a, t, x, y, v) by a disjunction describing all the sample 
points belonging to the trajectory sample a; the same for a'; and (2) eliminating all 
the quantifiers in the obtained formula. In concreto, using the notation from above, 
each occurrence of S{a,t,x,y,v) would be replaced in (y5aiibi[a, a'] by Vi!=o^(^ ~ 
ti Ax = Ay = YiAv ~ v^), and similar for a'. This results in a (rather complicated) 
first-order formula over the reals <^aiibi[a, a'] in which the predicate S does not occur 
any more. Since first-order logic over the reals admits the elimination of quantifiers 
(i.e., every formula can be equivalently expressed by a quantifier- free formula), we 
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can decide the truth value of (^aiibi[a, a'] by ehminating all quantifiers from this 
expression. In this case, we have to eliminate one block of existential quantifiers. 

We can however simplify the quantifier-elimination problem. It is easy to see, 
looking at (f) above, that ^(^aiibi[a, a'] is equivalent to 

V V ilJaUM [U , Xj , y,, , tj+1 , Xj+1 , y,+i , Vj , t'j , x^ , y'j , t^-^^ , x^_^i , y^^-^^^ , v^] , 

where the restricted alibi-query formula ■tpaHbiiU, Xi, yi, ti+i,Xi+i,yi+i,Vi, t'^,x'^,y'j, 
ij+i, a;^-|_]^, j/j^^i ''^'j) abbreviates the formula 

3t3x3y{U < t < t,+i A t;- < t < t^+i A (a; - x,)^ + {y - y,f < {t ~ U)^vf A 

(x - + (y - < {U+i - tfvf A 

(x - x'^f + {y- y]f <it- f^)\f A (x - x'^^.f + (y - y',+,f < {t'^+, tfvf) 

that expresses that two beads intersect. 
So, the instantiated formula 

''Pauu [tn X, , yi, t,+i , Xi+1 , y^+i , v., , t^ , x^ , y^ , t'^^-^ , x^^^ , y^^^ , v;^] 

expresses (f). To eliminate the existential block of quantifiers {3t3x3y) from this 
expression, existing software-packages for quantifier elimination, such as QEP- 
CAD [Hong 1990], Redlog [Sturm 2000] and Mathematica [Wolfram 2007] can be 
used. We experimented QEPCAD, Redlog and Mathematica to decide if several 
beads intersected. The latter two programs have a similar performance and they 
outperform QEPCAD. To give an idea of their performance, we give some re- 
sults with Mathematica: the computation of tpaiiui^, 0, 0, 1, 2, 2, VS, 0, 3, 3, 1, 2, 2, 
2] took 6 seconds; that of VaZ»bj[0, 0, 0, 1, 2, 2, a/8, 0, 3, 4, 1, 2, 2, 2] took 209 seconds 
and the computation of ipaiibi[^, 0, 0, 1,-1,-1, 1, 0, 1, 1, 2,-1, 1, 2] took 613 seconds. 
Roughly speaking, our experiments show that, using Mathematica , this quantifier 
elimination can be computed on average in about 2 minutes (running Windows 
XP Pro, SP2, with a Intel Pentium M, 1.73GHz, 1GB RAM). This means that 
evaluating the alibi query on the lifeline necklaces of two moving objects that each 
consist of 100 beads would take around 100 x 100 x 2 minutes, which is almost two 
weeks, when applied naively and at most (100 -I- 100) x 2 minutes or a quarter day, 
when first the intersection of time-intervals is tested. Clearly, in both cases, such 
an amount of time is unacceptable. 

There is a better solution, however, which we discuss next, that can decide if two 
beads intersect or not in a couple of milliseconds. 

3.4 The parametric alibi query 

The uninstantiated formula 

4'aHbi(ti, Xi, yi, Xi+i, J/i+l, Vi, tj, Xj, yj, ij+i, Xj^i, yj+l) Vj) 

can be viewed as a parametric version of the restricted alibi query, where the free 
variables are considered parameters. This formula contains three existential quan- 
tifiers and the existing software-packages for quantifier elimination could be used 
to obtain a quantifier- free formula ipaiibiiU, Xi, yi, ti+i,Xi+i,yi+i,Ui, t'^,x'j,yj,t'j_^_-^, 
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y'j+i^v'j) that is equivalent to ipaUbi- The formula ipaiiM could then be used to 
straightforwardly answer the alibi query in time linear in its size, which is indepen- 
dent of the size of the input and therefore constant. We have tried to eliminate the 
existential block of quantifiers 3t3x3y from tpaUM using Mathematica, Redlog and 
QEPCAD. After some minutes of running, Redlog invokes QEPCAD. After several 
days of running QEPCAD on the configuration described above, we have inter- 
rupted the computation without result. Also Mathematica ran into problems with- 
out giving an answer. It is clear that eliminating a block of three existential quan- 
tifiers from a formula in 17 variables is beyond the existing quantifier-elimination 
implementations. Also, the instantiation of several parameters to adequately cho- 
sen constant values does not help to produce a solution. For instance, without loss 
of generality we can locate {ti, Xi, yi) in the origin (0, 0, 0) and locate the other apex 
of the first bead above the j/-axis, i.e., we can take x^+i = 0. Furthermore, we can 
take Vi = 1 and t^+i — 1. But Mathematica, Redlog and QEPCAD cannot also not 
cope with this simplified situation. 

The main contribution of this paper is a the description of a quantifier-free for- 
mula equivalent to il)aiibi{ti.,Xi,yi,ti+i,Xi+i,yi+i,Ui,t'^,x'^,y'^,t'^j^^,x'^^^,y'^^^,v'j). 
The solution we give is not a quantifier-free first-order formula in a strict sense, 
since it contains root expressions, but it can be easily turned into a quantifier-free 
first-order formula of similar length. It answers the alibi query on the lifeline neck- 
laces of two moving objects that each consist of 100 beads in less than a minute. 
This description of this quantifier-free formula is the subject of the next section. 

4. PRELIMINARIES ON THE GEOMETRY OF BEADS 

Before, we can give an analytic solution to the alibi query and prove its correctness, 
we need to introduce some terminology concerning beads. 

4.1 Geometric components of beads 

Various geometric properties of beads have already been described [Egenhofer 2003; 
Kuijpers and Othman 2007; Miller 2005]. Here, we need some more definitions 
and notations to describe various components of a bead. These components are 
illustrated in Figure 2. In this section, let p = {tp,Xp,yp) and q = (tq,Xq,yq) be 
two time-space points, with tp < tq and let fmax be a positive real number. 

The bead B(p, g, Umax) is the intersection of two filled cones, given by the equa- 
tions {x ~ Xpf + {y~ ypY < {t- <p)^Wmax A tp <t and (x - Xq)'^ + (y - yq)^ < 
{tq — i)^Umax ^ t l^tq rcspcctively. The border of its bottom cone is the set of all 
points (t, X, y) that satisfy 

*C- (i, X, y, tp, Xp, yp, w^ax) := {x ~ Xpf + {y - ypf = {t - tpfv^^^^ A tp < t 

and is denoted by C~(p,Wniax) or C~ (ip, Xp, j/p, Umax); and the border of its upper 
cone is the set of all points (i, x, y) that satisfy 

'^C+{tiX,y,tq,Xq,yq,Vn^l,^) != (x - Xq)'^ + {y ~ yqf = [tq - O^^max ^ t < tq 

and is denoted by C+(g, Umax) or C+(tg, a;,, y,, Umax)- 

The set of the two apexes of B(p, g, Umax) is denotes rB(p, g, Umax), i-c, tB(p, g, 

Umax) = {P, q}- 
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We call the topological border of the bead B(p, Wmax) its mantel and denote 
it by 9B(p, Winax)- It can be easily verified that the mantel consists of the set of 
points (t, X, y) that satisfy 

\E'c' (t, X, ?/, tp^ Xp , yp, tq , Xq , yq , t'max) • tp t'^ tq /\ 

{2x{xp - Xq) +xl-xl + 2y{yp - yq) + yq ~ y^ < v'^^^2t{tp ~ tq) + t^ ~ t^h 
{x - Xp)^ + {y- yp)^ = {t- tp)^wLx V (a; - Xq)^ + (y - yq)^ = {tq - tfvl^^^ 

A 2x{xp - Xq) + a-2 - xl + 2y{yp - yq) + yl - y^ > {2t{tp - tq) +1]- i^)) . 

The first conjunction describes the lower half of the mantel and the second con- 
junction describes the upper half of the mantel. The upper and lower half of the 
mantel are separated by a plane. The intersection of this plane with the bead is an 
ellipse, and the border of this ellipse is what we will refer to as the rim of the bead. 
We denote the rim of the bead B(p, g, Wmax) by pB(p, q, Vmax) and remark that it is 
described by the formula 

^ p{t^ X, y, tp, Xp, T/p, tg, Xg, y^, ^^max) ■ 

(x - Xp)2 + (y ~ yp)2 = {t- tp)''vl^^ Atp<t<tqA 
2x(xp - Xq) +xl~xl + 2y{yp - yq) + y^ - yp = u^^x (2i(ip - tq) + - tp) . 

The plane in which the rim lies splits the bead into an upper-half bead and a 
bottom-half bead. The bottom-half bead is the set of all points (t, x, y) that satisfy 

^p; '^p: Up: ^q^ Xg, y^, 't^max) ■ 

(x - Xp)' + (y - yp)2 < (t - tp)2wLx A tp < t < A 
2x(xp - Xq) + x^ - x^ + 2y(yp - yq) + yl - yl < <ax (2i(ip - ig) + - ip) 

and is denoted by B"(tp, Xp, yp, tg, Xg, y,, Umax)- 

The upper bead is the set of all points (t, x, y) that satisfy 

^ (t ^ X ^ y ^ tp ^ Xp , yp, ig , Xg , yg , i^max) ■ 

(x - Xg)2 + (y - yg)2 < {tq - A tp < t < tg A 

2x(xp - Xg) + x^ - x^ + 2y(yp - yg) + y, - yp > u^^x (2i(ip - ig) + ^g - ip) 
and is denoted by B+(tp, Xp, yp, ig, Xg, yg, ■ymax)- 
4.2 The intersection of two cones 

Let (ti, Xl, yi, wi) and (t2, 2^2, y2, ^'2) be two bottom cones. A bottom cone, 
e.g., C~(ti, Xl, yi, wi), can be seen as a circle in 2-dimcnsional space (x, y)-space 
with center (xi, yi) and linearly growing radius (t — ti)ui as ti < t. 

Let us assume that the apex of neither of these cones is inside the other cone, 
i.e., (xi - xa)^ + (yi - ya)^ > {ti - t2)^vf V ii < t2 and (xi - xa)^ + (yi - y2)^ 
> (^1 — t2)'^V2 y t2 < ti. This assumption implies that at ti and t2 neither radius is 
larger than or equal to the distance between the two cone centers. So, at first the 
two circles are disjoint and after growing for some time they intersect in one point. 
Wc call the first (in time) time-space point where the two circles touch in a single 
point, and thus for which the sum of the two radii is equal to the distance between 
the two centers the initial contact of the two cones (ti, xi, yi, wi) and C^{t2,X2, 
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Fig. 2. A dissection of the bead Q{tp,Xp,yp,tq,Xq,yq,Vuis.yi)- 



y2,V2)- It is the unique point (t,x,y) that satisfies the formula 

^ic-{t,x,y,ti,xi,yi,vi,t2,X2,y2,V2) := ti < t At2 < t A 
{x - xif + {y- yif = {t- hfvl A {x - X2Y + {y - ya)' = {t - t2)'f| A 

((t - ti)vi + (t - t2)f2)' - (a^i - 2:2)2 + (2/1 - yz)'. 

The initial contact of two cones Q.'^ {ti,xi,yi,vi) and C"''(t2, X2, y2j ^'2) is given by 
the formula (t, x, y, ti, Xi, wi, ^2, 2^2, 2/2, ^^2) that we obtain from \E'/(7- by 

replacing in ti < tAi2 < i by f < ti At < ^2- We denote the singleton sets containing 
the initial contacts by IC(C~ (ii, xi, yi, wi), C~(t2, a;2, y2, ^2)) and IC(C+(ti, xi, yi, 
wi),C+(t2,a:2,y2,'y2))- 



{ti,xi,yi) 



initial contact — 

{t,x,y) 



V2{t~t2) 




Fig. 3. Intersecting cones and their initial contact (3-dimensionaI view on the left and 2- 
dimensional view on the right). 
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From the last equation in of the system in ^'/c- and ic+ , we easily obtain 



t = x2)~+{yi y2)+tivi+t2V2 ^ compute the other two coordinates {x,y) of 



the initial contact, we observe that for in the plane of this time value t, it is on the 
line segment bounded by (xi,yi) and (2:2,2/2) and that its distance from {xi,yi) 
is vi{t — ti) and its distance from (xi,yi) is W2(i — ^2)- We can conclude that the 
initial contact has {t, x, ?/)-coordinates given by the following system of equations 



This means that we can give more explicit descriptions to replace ^'/c- and 5'/c+. 
5. AN ANALYTIC SOLUTION TO THE ALIBI QUERY 

In this section, we first describe our solution to the alibi query on a geometric level. 
Next, we prove its correctness and transform it into an analytic solution and finally 
we show how to construct a quantifier-free first-order formula out of the analytic 
solution. 

5.1 Preliminary geometric considerations 



The solution we present is based on the observation that the two main cases of 
intersection (that do not exclude each other) are: (1) an apex of one bead is in the 
other; and (2) the mantels of the beads intersect. 

The inclusion of one bead in the other, illustrated in Figure 4. is an example of 
the first case. It is clear that if no apex is contained in another bead and we still 
assume that the beads intersect, than their mantels must intersect. We show this 
more formally in Lemma 1. In this second case, the idea is to find a special point 
(a witness point) that is easily computable and necessarily in the intersection. 





Fig. 4. One bead is contained in the other. 
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Fig. 5. Clean cut between cones. 

Let us consider two beads with bottom cones C~{ti,xi,yi,vi) and {t2, X2,y2i 
V2) and let us assume that none of the apeces is inside the other cone. One special 
point is the point of initial contact IC(C^ (ii, xi, wi), (<2, 2:2, 2/2, ^2)). However, 
this point can not be guaranteed to be in the intersection if the mantels of the two 
beads intersect, as we will show in the following example. Consider two beads with 
bottom cones C^(0, 0,0, 1) and C^(0, 2,0, 1). The intersection is a hyperbola in 
the plane x — 1 with equation — = 1. The initial contact of the two bottom 
cones is the point (1,0, 1). To show that this point of initial contact docs not need 
to be in the intersection of the two beads, the idea is to cut this point out of the 
intersection as follows. Suppose one bead has apexes, (0, 0, 0) and (a, 6, c) and speed 
1. The plane in which its rim lies is given by —2ax + — 2hy + 6^ + 2ct — = 0. 
This plane cuts the plane a given by the equality a; = 1 in a line given by the 
equation —2hy + 2ct — 2a + — =0. Clearly, we can choose (a, 6, c) such that 

the line contains the points 1, and (v^, 1, l)- Everything below this line 
will be part of the first bead and the second cone, but the initial contact is situated 
above the line, effectively cutting it out of the intersection. All this is illustrated 
in Figure 6. 

We notice how the plane in which the rim lies and the rim itself is the evil do- 
er. If neither rim intersects the mantel of the other bead, then the intersection of 
mantels is the same as an intersection of cones. In which case the initial contact 
will not be cut out and can be used to determine if there is intersection in this 
manner. 

Using contraposition on the statement in the previous paragraph we get: if there 
is an intersection and no initial contact is in the intersection then a rim must 
intersect the other bead's mantel. 

To verify intersection with the apexes and initial contacts is straightforward. 
Verifying if a rim intersects a mantel results in solving a quartic polynomial equation 
in one variable and verifying the solution in a single inequality in which no variable 
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Fig. 6. The initial contact cut out. 

appears with a degree higher than one. 
5.2 Outline of the solution 

Suppose, for the remainder of this section, we wish to verify if the beads Bi — 
B(ii,xi,yi,t2,X2,2;2,wi) and 82 = B(t3, 2:3, 2/3, i4, a;4, 1/4, ^2) intersect. Moreover, 
we assume the beads are non-empty, i.e., {x2 — a^i)^ + (2/2 — Ui)^ < (^2 ^ ti)^vl and 
(x4 - xs)'^ + (2/4 - ys)"^ < {ti - t3)'^vl 

We first observe that an intersection between beads can be classified into three, 
mutually exclusive, cases. The three cases then are: 

(I) an apex of one bead is contained in the other, i.e., 

tBi n B2 7^ or Bi n TB2 ^ 0; 

(II) not (I), but the rim of one bead intersects the mantel of the other, i.e., 

pBi n (9B2 7^ or pBj n dBi ^ 0; 

ACM Journal Name, Vol. V, No. N, 20YY. 



16 



(III) not (I) and not (II) and the initial contact of the upper or lower cones is in 
the intersection of the beads, i.e., 

IC(C]:, Q) c Bi n B2 or IC(Cj^, C^) c Bi n B2. 

If none of these three cases occur then the beads do not intersect, as we show in 
the correctness proof below. First, we give the following geometric lemma. 

Lemma 1. If Bi n B2 7^ 0, tBi n B2 = and rB2 n Bi = 0, then dBi n 882 ^ 0. 

Proof. From the assumptions, we know there is a point pi in B2, e.g., an apex 
of B2, that is not in Bi. Also, there is a point p2 that is in B2 and in Bi. The 
line segment bounded by pi and p2 lies in B2 , since B2 is convex. The line segment 
cuts the mantel of Bi since p2 is inside Bi and pi is not. Let p be this point where 
the segment bounded by pi and p2 intersects dBi. This point lies either on the 
upper-half bead B^ or on the bottom-half bead . Let r be the apex of this half 
bead. Since p is inside B2 and r is not, the line segment bounded by p and r must 
cut 9B2 in a point q. This point lies of course on dB2 and on dBi since the line 
segment bounded by p and r is a part of dBi. Hence their mantels must have a 
non-empty intersection if the beads have a non-empty intersection and neither bead 
contains the apexes of the other. □ 

Now, we show that if Bi and B2 intersect and neither (I), nor (II) occur, then 
(III) occurs. 

Theorem 1. If Bi n B2 7^ 0, tBi n B2 = 0, Bi n rB2 = 0, pB^ n ^62 = and 
pB2 n dBi = 0, then IC(Cf , C^) c Bi n B2 or IC(Cj^, C^) c Bi n B2. 

Proof. Let us assume that the hypotheses of the statement of the theorem is 
true. It is sufficient to prove that either nCj C B^ DB^ or C;^ nC^ C Bj^ n B^. 
We will split the proof in two cases. From the fourth and fifth hypotheses it follows 
that either (1) pB^ C B2 or pBj C Bi; or (2) pB^ n B2 = and pBa n Bi = 0. 

Case (1): We assume pBj C Bi (the case pB^ C B2 is completely analogous). We 
prove Cf n C Bj^ n B^ (the case for upper cones is completely analogous). The 
following argument is illustrated in Figure 7. 

Since pB2 C Bi, we know that pBj is inside Cj^, and (t^, X3, j/3) is outside. We can 
show that V2 < ui . Consider the plane spanned by the two axis of symmetry of both 
and C2 . Both and C2 intersect this plane in two half lines each. Moreover, 
we know that intersects the axis of symmetry of C^. Let Iq be the moment 
at which this happens. Obviously to > ti, but we know also know to > ts since 
{t3,X3,y3) is outside . We have that vi{to - ti) = ^ (xi - x^Y + (2/1 - 2/3)^- 
Since pBj is inside C]^ and (^3, 0:3, 2/3) is outside, this means both half lines from 
intersect the half lines from C]^ . Let t'^ and be the moments in time at which this 
happens and let t'^ > Iq. We have again that > ti and tj, > t3. Then vi{tQ — ti) = 
y/{xi - xs)^ + (yi - y3)^+W2(id - ^3) if and only if viit'^ - to) ^ V2{t'o - t^) . Since 
to > ^3, we get V2 < vi. This is depicted in Figure 7. 

It follows that every straight half line starting in [t^, X3, j/3) on C2 intersects Cj^ 
between (i3,a;3,j/3) and pB2, since pB2 is inside Cj^, and (t3, 0:3, 2/3) is outside. We 
also know that this line does not intersect beyond pBj since the cone is 
entirely inside Cj^ beyond the rim pBj. Therefore, n C2 C B2 ■ 
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Fig. 7. Illustration to the proof. 

Clearly, intersects since it can not intersect B^. We know Cf n OBj" is a 
closed continuous curve that lies entirely in . This curve is also contained in Bj^. 
Indeed, if we assume this is not the case, then it intersects the plane in which pBi 
lies, and hence it intersects pB^ itself, contradicting the assumption pB^ n dB2 ~ 0. 

Case (2): Now assume pB]^nB2 — and pEjOBi = 0. Clearly, vi cannot be equal 
to V2, otherwise the depicted intersection can not occur. So suppose without loss of 
generality that V2 < f i . Now either B^ intersects both B^ and B^ or B^ intersects 
both Bf and B^. These cases are mutually exclusive because of the following. If 
B^ intersects B^ then pBj is inside C^, likewise if B^ intersects B^ then pB2 is 
inside Cj^. Hence pBj C Bi which contradicts our hypothesis. If B^ intersects B-^ 
then pB2 must be outside and thus B2 must be as well, hence B^ intersects 
neither B^ nor B^ . Likewise, if B^ intersects B]^ then B^ can not intersect Bi. 

To prove that if Bj" intersects B^ then it also intersects B^ and if B2 intersects 
B^ then it also intersects B^ we proceed as follows (the case for B^ is analogous). 
Suppose B2 intersects B^, then B^nB^ C Bi, but pBj is outside Bi, that means B^ 
must intersect Bj'" since it can not intersect B^ anymore. This is the "what goes in 
must come out" -principle. Likewise, suppose B^ intersects Bj^, then Bj^nBj'" C Bi, 
but (t3,X3,y3) is outside Bi, that means B^ must intersect Bj^ since it can not 
intersect B^ anymore. 

So suppose now that B^ intersects both B-^ and Bj'" (the case for Bj" is com- 
pletely analogous). If B^ intersects B^ that means pBj is completely inside Cj^ 
and therefore that n C B^. We proceed like in the first case, we know 
that Cf n B^ is a closed continuous curve. This curve lies entirely in . If this 
curve is not entirely in B^ that means it intersects the plane in which pBi lies, 
and hence intersects pBi itself. But this is contradictory to the assumption that 
pB^ n dB2 = 0. □ 

In Theorem 1, we proved that if there is an intersection and neither rim cuts the 
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Fig. 8. Case (II) is not redundant: B(0, 0, 0, 2, 0, 2, 1.9) and B(0, 3, 0, 2, 3, 2, 1.9) seen from the top 
and the side. 

other bead's mantel and neither apex of a bead is contained in the other then there 
must be an initial contact in the intersection. Visualizing how beads intersect might 
tempt one to think there is always an initial contact in the intersection. There exist 
counterexamples in which there is an intersection and no initial contact is in that 
intersection. That means case (II) is not redundant. This situation is depicted in 
Figure 8. The beads are B(0, 0, 0, 2, 0, 2, 1.9) and B(0, 3, 0, 2, 3, 2, 1.9). 

It is clear that the initial contact of the bottom cones lies in the plane spanned 
by the axis of symmetry of those bottom cones, in this case this is the plane y ~ 0. 
The intersection of Figure 8 can be seen in Figure 9, where the two beads clearly 
have no intersection and thus no initial contact in the intersection. 

In the case of the upper cones the initial contact must lie in the plane y = 3. 
The intersection of Figure 8 can be seen in Figure 10, where the two beads clearly 
have no intersection and there is again no initial contact in the intersection. 
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Fig. 9. Intersection of Figure 8 with the plane y = 0. 




Fig. 10. Intersection of Figure 8 with the plane y = 3. 

This concludes the outline. 

5.3 A formula for Case (I) 

In Case (I), we verify whether rBi n B2 7^ or Bi n TB2 7^ 0. To check if that is 
the case we merely need to verify if one of the apexes satisfies the set of equations 
of the other bead. In this way wc obtain 

^"1 {tl,Xl,yi,t2,X2,y2,Vi,t3,X3,y3,t4,X4,y4,V2) : = 

(*B it3,X3,y3,ti,xi,yi,t2,X2,y2,vi) V *b (^4, 2:4, 2/4, ii, ^i, yi, ^2, 2^2, 2/2, I'l) V 

*B itl,Xi,yi,t3,X3,y3,t4,X4„y4,V2) V {t2,X2,y2,t3,X3,y3,t4,X4,y4:,V2)) ■ 

For the following sections we assume that the apex sets of the beads are not 
singletons, i.e., ti < t2 and ^3 < ^4. 

5.4 A formula for Case (II) 

Now, let us assume that $1 failed in the previous section. Note that we can 
always apply a speed-preserving [Kuijpers and Othman 2007] transformation to 
R X to obtain easier coordinates. We can always find a transformation such 
that {t'i,x[,y[) = (0,0,0) and that the line-segment connecting {t'i,x[,y[) and 
(^2, 0:2, 2/2) is perpendicular to the y-axis, i.e., y'2 — 0. This transformation is a 
composition of a translation in R x R^, a spatial rotation in R^ and a scaling in 
R X R^ [Kuijpers and Othman 2007]. Let the coordinates without a prime be the 
original set, and let coordinates with a prime be the image of the same coordinates 
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without a prime under this transformation. Note that we do not need to transform 
back because the query is invariant under such transformations [Kuijpers and 0th- 
man 2007]. The foUowing formula returns the transformed coordinates {t' , x',y') of 
{t,x,y) given the points (ti,xi,yi) and (^2, 2:2, j/2): 

ipA{ti,xi,yi,t2,X2,y2,t,x,y,t',x',y') ~ (y2 7^ J/i A 

t' = {t- ti)y/{x2 - xi)^ + {y2 - A x' = {x- xi){x2 - Xi) 

+iy - yi){y2 - yi) a y' ^ {x- xi){yi - 7/2) + (y - yi){x2 - xi)) 

V (2/2 = 2/1 A t'^{t-h) A = A = 

The translation is over the vector (— ti, — xi, — yi), the rotation over minus the 
angle that (t2 — ti^ X2 — xi,y2 — yi) makes with the x-axis, and a scaling by a factor 

{x2 — xi)'^ + (2/2 — 2/i)^- Notice that the rotation and scaling only need to occur 
if 2/2 is not already in place, i.e., if 2/2 7^ 2/i- 

The formula ilJcrdit[,x[,y[,ti,xi,yi,t'2,X2,y2,t2,X2,y2,t3,x'r^,y'^,t3,X3,y3,t'^,x'4, 
2/4, i4, 2:4, 2/4) is short for (fAlti,xi,yi,t2,X2,y2,ti,xi,yi,t[,x[,y[) A (pAih, xi,yi, 

t2,X2,y2,t2,X2,y2,t'2,X2,y'2) A (ti , Xi , 2/1 , ^2 , 2^2 , 2/2 , *3 , 2^3 , 2/3 , ^3 > 2:^3 > 2/3) A (fiAih, 
Xl,yi,t2, X2, 2/2, t4, X4, 2/4, ^4, ^^4, 2/4)- 

This transformation yields some simple equations for the rim pBi'. 



Not only that, but with these equations we can deduce a simple parametrization 
in the x-coordinate for the rim, 



We remark that this implies t2^Q and vi ^ 0. If = 0, then Bi is a point, hence 
degenerate. If vi = 0, then Bi is a line segment, and again degenerate. Next we 
will inject these parameterizations in the constraints for 9B^ and 9B^ separately. 
The constraints for 9B^ are 



{x ~ 4)2 + {y- y'jf = (i - t'^fvl 

2x{x'3 - x'i) + x'i - x'i + 2y{y'3 - 2/4) + 2/f - y'i < vl (2i(t(, - t'^) + t'i - t'i) 

t'3<t<t'^ . 



We will explain how to proceed to compute the intersection with ^i^d simply 
reuse formulas for intersection with 9Bj. First, we insert our expressions for x and 
2/ in the first equation. This is equivalent to computing intersections of pBj^ with 
and gives 





) 



2 




2^2 



) 



2 



-2:2-2/3 
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or equivalently 



-t' 



±2y',^Jvf {2xx', xq + viqf {IvH'^f x^ = 

(2x4 - ^'i + vft'i - {^v'A) t',)\i - {2vlt',f {x - x',r 

- {2vlt',yy'i - {vl {2xx', - x'i + vlt'if (2vlt',f x^) 
or equivalently 



±v,2y'^^xH [xq - vlt'i) + xAx'i [vlt'i ~ x'i) + {vlt'i - x'^f 
= xHx'i [vl - vl) + 
xA {^x'ivl {-x'q + vlt'i ~ Avtt'it':,) + 2vtt'ix':, + vlx'^ {vlt'i - x'i)) 
+ [vl {-xi + vlti - ivft'it'.Y ~ Avfti {x'i + y'i) - vt {^xi + vlt'i)) . 

By squaring left and right hand in this last expression, we rid ourselves of the 
square root and obtain the following polynomial equation of degree four. Squaring 
may create new solutions, so to ensure we only get useful solutions, we have to add 
the condition that the square root exists. This is the case if and only if 

cj)^{x, t'^,x'2, Vl) := xH {xi - vlt'i) + x4xi {vlti - x'i) + {vlti - x'if > 
is satisfied. 

We notice that if Bi is degenerate, i.e., x'i = vlt'i, then the square root vanishes 
and the polynomial in ^4 is the square of a polynomial of degree two, yielding to 
at most two roots and intersection points as we expect. The case were wi = is 
captured by the formula in the next section, that is why we leave that case out 
here and demand that vi ^ 0. So the following still works if one or both beads is 
degenerate: 

4>i{x, t'2, x'2, Vl, tg, Xg, j/g, 1)2) '■= 3a3b3c3d3e {ax"^ + bx^ + cx^ + dx + e ~ 

A a = {Axi {vi ~vl))' A b= -32x'ivl («| - vl) {^xi + vlti - Avft'it', 
+2vtt'ix', + vlx', {vlti - x'i)) A c = 8 {xi vlt'i) {-Av\t!i {x'i + y'i) + 
+vl {-xi + vlti - Avtt'it',)' - vt {~xi + vlt'i)) + {2viy',)' {xi - vlt'i) 

+ (4 {-x'ivl {-xi + vlti - Ht'it',) + 2vttix', + vlx', {vlti x'i)))' 
A d^8 {-xivl {-xi + vlti - Avft'it'^) + 2vttix'^ + vlx'^ {vlti - x'i)) 




ACM Journal Name, Vol. V, No. N, 20YY. 



22 



The quantifiers we introduced here are only in place for esthctical considerations 
and can be eliminated by direct substitution. 

We note that if vi = V2, we get polynomials of degree merely two. This can be 
solved in an exact manner using nested square roots (or Maple if you will). This 
gives us at most four values for x. Let 

4'roots{Xcn ^b: ^ci ^dj ^2 i -^21 ^1 1 ^3 i •''Si Vsi ^'2) 

be a formula that returns all four real roots, if they exist, that satisfy both (j^i^x, ^2' 
Vi,t'^, 2:3, 2/3, V2) and 4>^{x, X2, wi). We substitute these values in the parameter 
equations of pE^. By substituting these in the last equation above, we can deter- 
mine the sign of the square root we need to take for y. A point {t, x, y) satisfies the 
following formula is a point on pB]^, but instead of using the square root for y, we 
use an expression from above to get the correct sign for the square root if 2/3 7^ 0. 
If 7/3 = we have to use the square root expression and then it does not matter 
which sign the square root has; we need both: 

i'p{t,x,y,t'2,x'2,vi,t'^,x'-i,y'3,V2) := {y'3 ^ A t {2vlt'2) 2xx'2 ~ x'^ + I'lt'^ A 

2y's {2vlt'2) y = (2x4 " 4' + vftf - {^^2) t'zf vl - {^^2^ {x - 4)' 
- {2vlt',fy'i - (vl (2x4 - X2 + vlt'if - {2vlt',fx') A Q<t<t',) 
V (yg = A t {2vlt2) = 2a;a;2 - x'^ + vlt'^ A < t < A 
{2vlt'2f y" = (2x4 - 4' + vlt'i)^ - (2^'?4)' 4) . 

The four roots give us four spatio-temporal points on H Q.^- In order for 
these points (t, x, y) to be in pBj^ n 9B2', they need to satisfy 

V'-(t,x,y,t3,X3,y3,i4,X4,y4,W2) := 

2x(x^ - x^) + x'l - x'i + 2y(y^ - y^) + yf - y'i < i| (2i(4 - t',) + tf - t',^) . 

This formula returns True if (t, x, y) lies in the same half space as the bottom-half 
bead. 




Fig. 11. The rim intersects the eone and solutions arc verified in a half-space. 

The formula returns True if {t, x, y) lies in the same half space as the 
upper-half bead, i.e., ^+{t,x,y,t'^,x'^,y'r^,t'^,x'^,y'^,V2) ip-{t,x,y,t'^,x'^,y'4^,fr^,x'^, 
y3^V2). By combining ipp{t,x,y,t'2,X2,Vi,t'^,X3,y3,V2) and il;^{t, x,y,i, x,y,i, x,y, 
v) we get a formula that decides the emptyness of the intersection pB^ D in 
terms of a parameter x: 

V'pnai (2;, ^2, X2, wi, ^3, X3, yg, V2, i, x, y, i, x, y, v) := 
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3y (y^ = A {2vlt',f = {2xx'2 - x'i + vlt'^fvl - {2vlt'^f x^) V 
(y^ ^ A 2y'^ {2vlt'^) y - (2x4 ^ 4' + «242 _ (2«?4) t^)' - (2z;2t'2)' 

{X - - (2«2,. )2 ^,2 _ ^^2 (2.^,^, _ ^,2 ^ ^2,,2)2 _ (2^2,Q^2^^ ^ 

{2x{x - £) + ^2 _ f 2 ^ 2y(y - y) + - {2vlt'^) < 

(2 (2uji9 {2xx2 - xf + vjt'i) {i - i) + (2i;?iQ {P -~ P)) A < t'^ 

(2x4 - ^2 + vlt'^) < 2vlt2 A (£ (2vitf ) < 4 (2a;4 - 4^ + ^i^'i) < * (2w?if ) 

V i(2vlt'2) < t'2 {2xx'2 - 4^ + f'l^^) < H'^vft'i)) . 

We are ready now to construct the formula that decides if pBi and B2 have a 
non-empty intersection: 

'Pp,nd-(t'2^^'2iVi:t3:x'-i,y'2,t'4^,x'^,y'^,V2) -.^ 3x3xa3xb3xc3xd ( 

(l)rootsiXa,Xb,Xc,Xd,t'2,x'2,Vi,t'^,x'3,y'^,V2) A (x = V X = Xfc V X = Xc V X = .Trf) 
A V'pna±(2^,*2j47'f'l>*Js'4'2^3'*^2,i3,4'y3'*47 2^47y47'f'2) ) • 

The formula that decides if pBi intersects looks strikingly similar: 

'Ppind+(^'2'^'2''"iit'^^x'^,y'r^,t'4^,x'^,y'^,V2) -.^ 3x3xa3xb3xc3xd ( 

4>roots{Xa,Xb,Xc,Xd.,t'2.,x'2,Vi,t'^,x'4^,y'^,V2) A (x = Xa V X = X(, V X = Xc V X = Xd) 

A '(/'pna±(a;,t2,X2,'yi,t4,X4,y4,W2,t4,X4,y4,t3,X3,?;3,4) ) • 

The quantifiers introduced here can also be eliminated in a straightforward man- 
ner. Notice that (proots acts as a function rather than a formula that inputs 
(^2, 2^2; ''^i; "^2) to Construct a polynomial of degree four and returns the 

four roots (Xa, Xb, Xc, Xd), if they exist, of that polynomial. The existential quan- 
tifier for the variable x is used to cycle through those roots to see if any of them 
does the trick. Finally we are ready to present the formula for Case (II) : 

*ii (^I,a:i,yi,t2,a;2,y2,wi,i3,a;3,y3,i4,a:4,y4,t'2) := (wi 7^ A W2 7^ 0) A 
$1 (ti,xi,yi,t2,2:2,y2,wi,<3,a;3,y3,<4,a;4,y4,W2) A 
3t\3x\3y[3t'23x23y'23t'^3x'^3y'^3t'^3x'43y'^ ( 
■il^crdit[,x[,y[,ti,xi,yi,t'2,x'2,y'2,t2,X2,y2,t'3,x'ri,y'3,t3,X3,y3,t'^,x'^,y'4^,t4,X4,y4) 

^ (^</'pina-(*2:a;2>^^i:*3:a;3,y3,t4>a:4,y4,f2) V 

V'pina+ (*27 4: "1, *3> 4' y3> t'i, x\, y'^, V2)j V 
i^crd{t3,x'3,y3,t3,x3,y3,t'^,x'i,y'i,t4,x4,y4,t[,x[,y[,ti,xi,yi,t2,x'2,y'2,t2,x2,y2) 

^ ['Pp,nd-ih^X3,V2,h,X4,y4,t2,X2,y2,V4) V 

'^Pi na+ (^3 ' 2^3 : f 2 , i'l , a;'i , y'l , , 2^2 > y2 > «i )) ) • 

The reader may notice that a lot of quantifiers have been introduced in the for- 
mula above. These quantifiers are merely there to introduce easier coordinates and 
can be straightforwardly computed (and eliminated) by the formula ipcrd and hence 
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the formula (pA{ti.,xi,yi,t2,X2,y2,t,x,y,t' ,x' ,y'). The latter actually acts like a 
function, parameterized by (ii, xi, yi, t2, 2^2, 2/2), that inputs (t,x,y) and outputs 
{t\x',y'). 

5.5 A formula for Case (III) 

Here, we assume that both ipi and tpu fail. So, there is no apex contained in the 
other bead and neither rim cuts the mantel of the other bead. 

As we proved in Theorem 1, the intersection between two half beads will reduce 
to the intersection between two cones and that means there is an initial contact 
that is part of the intersection. To verify if this is the case we compute the two 
initial contacts and verify if they are effectively part of the intersection. 

Using the expression for the initial contact IC(C]^,C^), we computed in Sec- 
tion 4.2 we can construct a formula that decides if it is part of n B^. We 
will recycle the formulas V'- from the previous section to construct an expression 
without the need for extra variables. The following formula that returns True if 
IC(Cf ,C^) = {to,xo,yo) satisfies tl,'^{to,xo,yo,t' ,x' ,y' ,i,x,y,v): 

<j>- {ti,xi,yi,vi,t3,X3,y3, V2,t\ x' , y', i, x, y, v) := 

2{x' - x) (^{xiV2 + X3Vi)y^{xi - X3)2 + {y^ - y-^f + Vi {{t^ - ti)v2) {x^ - Xi)^ 

+ 2(y' - y) [{yiV2 + yzvi)\/{xi - x-iY + {yi - ys)^ + vi {{t^ - ti)v2) [ys - yi)j 
+ V{^i-^3Y + {yi-y3yi^i+v2) {x^ - x'^ + f- y'^) < {{t' - 1'^) («i + V2) 
+ 2 (^^{xi - xsY + (2/1 - 2/3)2 + + ^3^2) (t'-t)) Vi^i - + (yi - 2/3)' • 

The following formula expresses that the time coordinate to of IC(C|f, C^) satisfies 
the constraints t' < t < t" and t < t < i: 

ipt {ti,xi,yi,vi,t3,X3,y3,V2,t',t" ,i,t) := 

t\vi + V2) < \J {xx - X3Y + (2/1 - 2/3)2 ^ i^y^ ^ ^3^2 < i!' (v\ + V2) 
A i[vx + V2) < V (xi - xsY + {yi - 2/3)2 _|_ ^^y^ ^ ^3^2 < i{vi + V2) ■ 

Now, IC(Cj;,C2") C Bj^ n if and only if V'lC- (^ii 2:1, 2/1, ^2, a;2, 2/2, f 1, ^3, a^3, 2/3, 
^4, 2^4, 2/4,^2) where V'lC- (^i, a^i, 2/i, *2, 2:2, 2/2, wi, ^3, 2:3, 2/3, ^4, ^4, 2/4, ^2) := i>t{ti,Xi, 
yi,vi,t3,X3,y3,V2,ti,t2,t3,t4) A (f>^{ti, xi, yi, Vi,t3, X3, y3, V2, h, Xi, yi, t2, X2, y2, 
vi) A xi, 2/1, ■yi,t3,a;3, 2/3, W2,t3,a;3, 2/3, t4,a;4, 2/4,^2) and iC(C^,C^) C B]|" n 

B^ if and only if ^ic+ (^i, a;i, 2/1, ^2, ^2, 2/2, t"!, ^3, 2:3, 2/3, ^4, 2:4, 2/4, 1^2) := A{t2,X2, 

2/2,^1,^4, 3^4, 2/4, 1^2, tl,i2,i3,i4) A (^2 , 2^2 , 2/2 , ""l , ^4, 2^4 , 2/4, ^2 , ^2 , 2^2 , 2/2 , , , 2/1 , 

Vi) A 0_(t2, 2:2, 2/2, fl,i4, 2^4, 2/4,^2,^4, 2^4, 2/4,^3, 2:3, 2/3, ^'2). 

The formula that expresses the criterium for Case (III) then looks as follows: 

$111 itl,Xi,yi,t2,X2,y2,Vl,t3,X3,y3,t4,X4,y4:,V2) ("l + V2 ^ 0) A 
^ $1 (tl, 2:1, J/1, t2, 2:2, 2/2,^1,^3, 2:3, 2/3,^4, 2:4, 2/4,^2) A 
(V'lC- 2^1, 2/1, ^2, 2:2, 2/2, t^l, ^3, 2^3, 2/3, h, X4, 2/4, V2) V 
V'|C+ (il, 2^1, 2/1,^2, 2:2, 2/2, t'1,^3, 2:3, 2/3,^4, 2:4, 2/4,^2)) . 
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5.6 The formula for the parametric alibi query 

The final formula that decides if two beads, Bi = B{ti, xi,yi,t2,X2,y2,vi) and 
B2 = B(i3, xs, ?/3, t4, X4, ?/4, W2), do not intersect looks as follows 

i^aiibi{ti,Xi,yi,t2,X2,y2,vi,h,X3,y3,t4,X4,y4,V2) --^ ^{{ti <t2 A ^3 < ^4) A 

($111 {tl,Xi,yi,t2,X2,y2,Vl,t3,X3,y3,t4,X4,y4,V2) 

V $11 {ti,Xi,yi,t2,X2,y2,Vi,t3,X3,y3,t4,X4,y4,V2)) 

V $1 {ti,Xi,yi,t2,X2,y2,Vl,t3,X3,y3,t4,X4,y4,V2)) ■ 

6. EXPERIMENTS 

In this section, we compare our solutionto the alibi query (using the formula given 
in Section 5.6) with the method of eliminating quantifiers of Mathematica. 

In the following table it is clear that traditional quantifier elimination performs 
badly on the example beads. Its running times highly deviates from their average 
and range in the minutes. Whereas the method described in this paper performs 
in running times that consistently only needs milliseconds or less. This shows our 
method is efficient and our claim, that it runs in milliseconds or less, holds. 

For this first set of beads we chose to verify intersection of two oblique beads 
(1-2) and the intersection of one obhque and one straight bead (3-4). The beads 
that actually intersected had a remarkable low running time with the QE-method. 





The beads 


The running times 


Bi 


B2 


QE 


Our Method 


1 

2 
3 

4 


(0,0,0,2,0,2,1.9) 
(0,0,0,2,0,2,1.9) 
(0,0,0,2,0,2,1.9) 
(0,0,0,2,0,2,1.9) 


(0,3,0,2,3,2,2) 
(0,4,0,2,4,2,2) 
(0,3,0,2,3,0,2) 
(0,4,0,2,4,0,2) 


0.656 Seconds 
324.453 Seconds 

0.438 Seconds 
475.719 Seconds 


0.016 Seconds 
0.063 Seconds 
0.015 Seconds 
0.031 Seconds 



The type of beads in this second set are as in the first. However, these beads 
all have overlapping time intervals unlike the first set, where the time intervals 
coincided. 





The beads 


The running times 




Bi 


B2 


QE 


Our Method 


1 

2 

3 

4 


(0,0,0,2,0,2,1.9) 
(0,0,0,2,0,2,1.9) 
(0,0,0,2,0,2,1.9) 
(0,0,0,2,0,2,1.9) 


(1,3,0,3,3,2,2) 
(1,4,0,3,4,2,2) 
(1,3,0,3,3,0,2) 
(1,4,0,3,4,0,2) 


63.375 Seconds 
59.485 Seconds 
29.734 Seconds 
27.281 Seconds 


0.078 Seconds 
0.078 Seconds 
0.031 Seconds 
0.032 Seconds 



The type of beads in this third set are as in the first. But this time the time 
intervals are completely disjoint. Note that the running times for the QE-mcthod 
are more consistent in this set and the previous one. 





The beads 


The running times 


Bi 


B2 


QE 


Our Method 


1 


(0,0,0,2,0,2,1.9) 


(3,3,0,4,3,2,2) 


63.641 Seconds 


0.046 Seconds 


2 


(0,0,0,2,0,2,1.9) 


(3,4,0,4,4,2,2) 


61.781 Seconds 


0.016 Seconds 


3 


(0,0,0,2,0,2,1.9) 


(3,3,0,4,3,0,2) 


52.735 Seconds 


0.031 Seconds 


4 


(0,0,0,2,0,2,1.9) 


(3,4,0,4,4,0,2) 


56.875 Seconds 


0.046 Seconds 
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7. THE ALIBI QUERY AT A FIXED MOMENT IN TIME 

7.1 Introduction 

In this section, we present another example where common sense prevails over the 
general quantifier-elimination methods. The problem is the following. As in the 
previous setting, we have lists of time stamped-locations of two moving objects and 
upper bounds on the object's speed between time stamps. We wish to know if two 
objects could have met at a given moment in time. 

For the remainder of this section, we reuse the assumptions from the previous 
section. We wish to verify if the beads Bi = B(ti, xi, yi, ^2, 2^2, ?/2, f i) and B2 = 
B(i3, X3, ?/3, t4, X4, ^4, W2) intersect at a moment in time to. Moreover, we assume 
the beads are non-empty, i.e., (^2 — xi)^ -f (?/2 — J/i)^ < {t2 — ti)^vfa.nd {x^ — xs)"^ + 
(2/4 — 2/3)^ < (*4 — HYv2 and that < to < ^2 and t3 < to < ^4 are satisfied. This 
means we need to eliminate the quantifiers in 

3x3y {{x - xif + {y- yif < vfito - hf A {x ~ X2)' + {y - 1/2)' < w?(to - ^2)' 
Aix- xsf + {y~ ysf < vl{to - tsf A {x ~ x^f + {y - y^f < vl{ta - tif) . 

Eliminating quantifiers gives us a formula that decide whether or not four discs 
have a non-empty intersection. For ease of notation we will use the following ab- 
breviations: € Di if and only if [x — XiY + {y — yi) < and {x,y) G Ci if 
and only if {x - x^)"^ + {y - Vi) = rf. 

7.2 Main theorem 

Using Helly's theorem we can simplify the problem even more. Helly's theorem 
states that if you have a set S oi m convex sets in n dimensional space and if any 
subset of of n 4- 1 convex sets has a non-empty intersection, then all m convex sets 
have a non-empty intersection. For the plane, this means we only need to find a 
quantifier free-formula that decides if three discs have a non-empty intersection. For 
the remainder of this section assume that we want to verify whether Di n D2 H D3 
is non-empty. 

Theorem 2. Three discs, Di, D2 and D3, have a non-empty intersection if and 
only if one of the following cases occur: 

(1) there is a disc whose center is in the other two discs; or 

(2) the previous case does not occur and there exists a pair of discs for which one 
of both intersection points of their bordering circles lies in the remaining disc. 

Proof. The ^/-direction is trivial. The only i/-direction is less trivial. We will 
use the following abbreviations, D = Di n £'2 H 1)3 and C = dD. 

Assume D is non-empty and that neither (1) nor (2) holds. The intersection D 
is convex as it is the intersection of convex sets. We distinguish between the case 
where Z? is a point or and the case where D is not a point. 

— Suppose D is a single point p. This point p can not lie in the interior of the three 
discs, because D would not be a point then. 

Nor can p lie in the interior of two discs. If that would be the case then there 
exists a neighborhood of p that is part of the intersection of those two discs, 
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say Di and D2 ■ Moreover p would be part of C3 and this neighborhood would 
intersect the interior of D3. This means D is not a point. 

So p must lie on the border of two discs, say Di and D2, and p must also be part 
of -D3 because D = {p}. This contradicts our assumption that (2) does not hold. 

— Assume D is not a point. All points on C belong to at least one C,;. If there is a 
point that docs not belong to any Ci , then it is in the interior of all Di and there 
exists a neighborhood of that point that is in the interior of all Di and hence in 
D. That contradicts to the fact that this point is in C. 

Furthermore, not all points of C belong to a single Ci. If that was the case then 
Di would be part of (and equal to) D and its center would be inside the other 
two discs which contradicts the assumption that (1) does not hold. 
So, C is made up of parts of the d, of which some may coincide but not all of 
them. When traveling along C you will encounter a point p that connects part of 
a Ci and part of a Cj, where i ^ j, that do not coincide, otherwise (1) must occur 
again which is a contradiction. However, this p also yields to a contradiction since 
it belongs to two different Ci, say Ci and C2, and is part of C hence D and D^. 
This contradicts the assumption that (2) does not occur. □ 

7.3 Translating the theorem in a formula 

We can simplify the equations even further using coordinate transformations. By 
applying a translation, rotation and scaling we may assume that {xi,yi) = (0,0), 
2^2 ^ 0; ri = 1 and 1/2 — 0. Using these simplifications and translating Theorem 2, 
we get the following formula. 

'fiix2,r2,X3,y3,r3) := {{-X2f <rl A (-X3)' + (-yg)' < V 

3x3y{x^ + i/ = l A {x - X2 f + ri A {x - xs)^ + {y - y^ f < r^) . 

This is a formula that decides if either the center of the first disc is part of the two 
other discs, sec the first line, or if either there exists a point in the intersection of 
the first two circles that is part of the third disc. All that remains now is making 
the expression 

3x3y{x^ + y^^l A {x ~ X2f + y^ ^ A {x - x^)^ + {y ^ y^f < rl) 
quantifier free. 

To do this we assume that Ci and C2 do not coincide but have a non-empty 
intersection. This is equivalent to 2:2 ^ A a;2 < ^2 + 1. Next, wc need to compute 
the point (s) where Ci and C2 intersect. 
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If y = 0; then verifying if that single point of intersection is part of D3 is easy, one 
only needs to verify if 

or equivalently (l — f2 ~ ^2)^ + ^^21/3 ^ '^^2^3 ■ 
If y 7^ 0, then verifying if one both points of interseetion is part of D3 is less trivial, 
sinee this involves square roots 

(^^^ - + (*2^^('-'-->'")((' + -)^--")- '■^'<- 
{xl + l~ rl ~ 2x2x3)'+ |^±y {rl - {X2 - if) ((1 + X2f - ri) - 2x22/3 j < ^xlrl 

or (x^ + 1 - - 2x2x3)' + {rl - (x2 - 1)^) ((1 + xa)^ - r^) 
+ (2x2y3)' ± 4x2^3^ (rl - (x^ - if) ((1 + X2)' - r^) < ixlrj 

or (x^ + 1 - - 2x2x3)' + (rl - (x2 - 1)') ((1 + X2)' - rl) 

+ (2x2^3)^ - Axlrj < ±4x22/3 y - (x2 - 1)') ((1 + X2)' - r|) . 

This is almost a F0(+, x, <, 0, l)-forinula except for the square root. However, 
the square root can be eliminated as we will show next. The previous expression 
is of the form L < ±a\/W. The presence of the ± simplifies this a lot, this means 
either sign of the square root will do, and also that we may assume the right hand- 
side is positive. Of course the square root must exist as well, this means > 0. 

This expression can then be simplified to 

W>0A{L<0\/L^< a^W) 
and gives us the expression 

$2(^2, r2,X3, 2/3,^3) (^2 - (2:2 - 1)') [il + X2f-rl) >0 

A ((x2 + l-rl~ 2x2x3)' + {rl - (X2 - 1)') ((1 + X2)' - rl) + (2x22/3)' 
- Axjrl < V ((x^ + l-rl- 2x2X3)' + {rl - (x2 - 1)') ((1 + X2)' - rf) 
+ (2x22/3)' - ixlrlf < (4x22/3)' {rl - {x2 if) ((1 + ^2)' - rl)) . 
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7.4 The safety formula 

Now, all that remains to be constructed is a formula that returns the convenient 
coordinates and a formula that guarantees that Ci and C2 actually intersect for 
safety, i.e., to exclude the case of empty intersection. The latter is constructed as 
follows. The formula f{xi,yi,ri,X2,y2,'f2) returns True if and only if the two 
circles, with centers {xi,yi) and {x2,y2) and radii ri and r2 respectively, have a 
distance between their centers that is not larger that the sum of their radii and not 
equal to zero to ensure they do not coincide. We have 

ip{xi,yi,ri,X2,y2,r2) := < {x2 - xif + {y2 ~ yif < [n + r2f ■ 

The formula 4>{xi,yi,ri,X2,y2,T2) returns True if and only if the second circle is 
not fully enclosed by the first, i.e., the sum of the distance between the centers plus 
the second radius is bigger than the first radius and vice versa. We can write 

4>{xi,yi,ri,X2,y2,r2) := {x2 - xif + (?/2 - yiY > ('^i - ''2)^ • 

These two safety conditions give us our safety formula 

*safety(^i'2^i'''i'^2,2;2,r2) := (p(xi,?/i,ri,X2,?/2,r2) A (l){xi,yi,ri,X2,y2,r2) . 

7.5 The change of coordinates 

The transformation consists of a translation, rotation and scaling. The translation 
to move the first circle's center to the origin. The rotation to align the second 
center with the x-axis. Finally the scaling to ensure that the first circle's radius is 
equal to one. First, the translation T{x, y) := {x — xi,y — yi). The rotation is 



R{x,y) 



^{x2~xiY + {y2-yiY \y1-y2 X2-X1J \ y 



X2~ xi y2- yi \ X 



and finally, the scaling is S{x,y) := y). The transformation is then a compo- 
sition of those three transformations A{x, y) = [S o Ro T)(x, y) := 

( {x2 - xi){x - xi) + {y2 - yi){y - yi) [yi - y2){x - xi) + {x2 - xi){y - yi) \ 
\ riyj{x2 - xiY + (7/2 - yif ' riyj{x2 - xiY + (j/2 - yiY j 

The following formula takes three circles with centers (xi,yi) and radii respec- 
tively, and transforms them in three new circles where the first circle has center 
(0, 0) and radius 1, the second circle has center {x2, 0) and radius r2 and the third 
circle has center {x'^,y'^) and radius rg. 

*transformation(2^i , 2/i , J"! , 2^2 , 2/2 , ^2 , ^3 , 2/3 , ^-3 , 2:2 , ?-2 , 2:3 > 2^3 : '"s ) : = 

^, ^ V{x2-x,r + iy2- ml r_2 ^ r,^ 

n ' n ^ n 

I {x2 - xi){xs - xi) + {y2 - yi){y3 - yi) , 



V3 = 



r\yj(x2 -xif ^ (2/2 - y\f 
[yi - y2){x3 - xi) + {x2 - xi){y3 - yi) 
ri\/{x2 - xiY + {y2 - yiY 
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Note that this is not a F0(+, x, <,0, l)-formula anymore due to the square roots 
and fractions. This "formula" is meant to act hke a function, which substitutes coor- 
dinates. The substituted coordinates have fractions and square roots but these can 
easily be disposed of when having the entire inequality on a common denominator, 
isolating the square root and squaring the inequality, as we showed in Section 7.3. 

7.6 The formula for the alibi query at a fixed moment in time 

First, we construct a formula that checks for any of two circles out of three if any 
of the conditions in Theorem 2 are satisfied. 

*2/3(^i'yi'''i'^2,y2,?'2,a;3,y3,''3) := ^x'r^Br'^Bx'^By'r^Br'r^ { 

^transformation (^^i > . : ^2 , 2/2 , r2 , X3 , ys , , , , , 2/3 , ) A 

(^*i(4,r^,4,?;^,r^) V $gg^fg|.y(xi, yi, ri, xs, ?/2, A $2(4. ^2. 4' ^3)) 

V *transformation(^i'yi'''i'^3,y3,7-3,a;2,y2,r2,a;2>''2:a;3,?/3,r^) A 
{^iix'2,r'2,x'3,yii,r'^) V $gg^fg|-y(a;i, ri, 0:3, 2/3, ra) A $2(2:2, 7-2, a^3' ^3)) 

V ^transformation ('^2 , ^2 , ''2 , a^s , ^3 , '^3 , a^i , yi , ''i , a;2 , '"2 , 2:3 , 2/3 , ) A 

(^5'i(4,r^,4,?/^,r^) V «'gafcty(2:2,2/2,?'2,a:3,2;3,f3) A «'2(4> ^2, 2:3, 2/3- ''3)) ) 

This formula is all we need to incorporate Helly's theorem in our final formula. 
Four discs have a non-empty intersection if and only if the following formula is 
satisfied 

^'(2:1, 2/1, J-i, a;2, 2/2, r2,X3,y3, rs, a;4, 2/4, r^) := 

*2/3(^l'yi'^l'^2,2/2,?-2,2;3,2/3,7-3) A «'2/3(2:i,2/l,ri,X2,2/2,7-2,2:4,2/4,?'4) 
A *2/3(^l' yi'''!' ^3, 2/3, ?'3, 2:4, 2/4,^4) A *2/3(^2,2/2,7'2,2;3,2/3,7'3,2;4,2/4,?'4)- 

This is almost a quantifier free-formula except for the fractions and square roots. 
However, as we showed before these can easily be disposed of. We omitted these 
tedious conversions for the sake of clarity. 

8. CONCLUSION 

In this paper, we proposed a method that decides if two beads have a non-empty 
intersection or not. Existing quantifier-elimination methods could achieve this al- 
ready through means of quantifier elimination though not in a reasonable amount 
of time. Deciding intersection of concrete beads took of the order of minutes, while 
the parametric case could be measured at least in days if a solution would ever 
be obtained. The parametric solution we laid out in this paper only takes a few 
milliseconds or less. 

The solution we present is a first-order formula containing square root-expressions. 
These can easily be disposed of using repeated squarings and adding extra condi- 
tions, thus obtaining a true quantifier free-expression for the alibi query. 

We also give a solution to the alibi query at a fixed moment in time. 

The solutions we propose are based on geometric argumentation and they illus- 
trate the fact that some practical problems require creative solutions, where at least 
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in theory, existing systems could provide a solution. 
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Appendix: The Mathematica implementation 

The function alibi, using the method described in this paper, returns True if the 
bead with apexes (tl, x\, yl) and {t2, x2, y2) and speed vl intersects the bead with 
apexes {t3,x3,y3) and (M,a;4,?/4) and speed v2 and False otherwise. The function 
alibiQE does the same except it uses the buih-in quantifier ehmination method of 
Mathematica. 

Note that Case (II) in our implementation corresponds to Case (III) in the 
description and vice versa. The reason for doing so is that Case (III), in the 
description, is computationally a lot easier than Case (II). Moreover, once any of 
the three cases returns True, our implementation exits and returns that result and 
thus omitting further computation. 



ACM Journal Name, Vol V, No. N, 20YY. 



33 



Needs [ "Miscellaneous" RealOnly" " ] ; 
isBeadNonEmpty[tl_, xl_, yl_, t2_, x2_, y2_, v_] : = 
tl < t2 A (x2 -xl) *2 + {y2 -yl) *2 < v'2* (t2 -tl) "2 
inbead[t_, x_, y_, tl_, xl_, yl_, t2_, x2_, y2_, v_] : = 
(x-xl) ^2 + (y-yl) *2 < v*2* (t - tl) *2A 
(x-x2) *2 + (y-y2) ^2 < v*2* (t -t2) *2 A (tl < t < t2) 
Casel[tl_, xl_, yl_, t2_, x2_, y2_, vl_, t3_, x3_, y3_, t4_, x4_, y4_, v2_] : = 
inbead[t3, x3, y3, tl, xl, yl, t2, x2, y2, vl] V 
inbead[t4, x4, y4, tl, xl, yl, t2, x2, y2, vl] Vinbead[tl, xl, yl, t3, 
x3, y3, t4, x4, y4, v2] Vinbead[t2, x2, y2, t3, x3, y3, t4, x4, y4, v2] 
CheckICCondition[tl_, xl_, yl_, vl_, t3_, x3_, y3_, v2_, lc_] : = 
- ({x3-xl) *2+ {y3-yl) ^2 < vl^2* {t3-tl) *2 Alf [Ic, t3 i tl, t3 s tl]) A 

- {(x3-xl) ^2 + (y3-yl) ^2 S v2*2* (t3-tl) -^2 Alf [Ic, t3 < tl, t3 > tl]) A 

- (vl == Av2 == 0) 

CheckHalf Space [t_, x_, y_, tl_, xl_, yl_, t2_, x2_, y2_, vl_, lc_] : = 
If [Ic, (tl < t S t2) A2*x* (xl -x2) +x2*2 -xl^2 + 2*y* (yl -y2) +y2*2 -yl^2 < 
vl'^2* (2*t* (tl-t2) +t2^2-tl^2), 
(tl <t <t2) A2*x* (xl-x2) +x2*2-xl*2 + 2*y* (yl-y2) +y2^2-yl*2 > 
vl*2* (2*t* (tl-t2) +t2*2-tl^2)] 
Case2[tl_, xl_, yl_, t2_, x2_, y2_, vl_, t3_, x3_, y3_, t4_, x4_, y4_, v2_] : = 
Module [{Ics}, If [CheckICCondition[tl, xl, yl, vl, t3, x3, y3, v2. True], 

Module [{to, xO, yO), tO = (tl * vl + 13 * v2 + Sqrt [ (xl - x3) ^2 + (yl - y3) ^2] ) / (vl + v2) ; 
xO = xl + vl* (tO-tl) * (x3-xl) /Sqrt[(xl-x3) ^2+ (yl - y3) *2]; 
yO = yl + vl* (tO-tl) * (y3-yl) / Sqrt [ (xl - x3) ^2+ (yl-y3) *2]; 
Ics = CheckHalf Space [to, xO, yO, tl, xl, yl, t2, x2, y2, vl. True] A 
CheckHalf Space [to, xO, yO, t3, x3, y3, t4, x4, y4, v2. True]], Ics = False] ; 
If [Ics, Ics = True, If [CheckICCondition[t2, x2, y2, vl, t4, x4, y4, v2. False], 
Module [{to, xO, yO}, tO = (t2*vl + t4*v2 + Sqrt[(x2 -x4) *2 t (y2 - y4) *2] ) / (vl + v2) ; 
xO = x2 + vl* (to -t2) * (x4 -x2) / Sqrt [ (x2 - x4) *2 + (y2 -y4) ^2] ; 
yO = y2 + vl* (t0-t2) * (y4-y2) / Sqrt [ (x2 - x4) ^2+ (y2-y4) ^2]; 
Ics = CheckHalfSpace[tO, xO, yO, tl, xl, yl, t2, x2, y2, vl. False] A 

CheckHalfSpace[tO, xO, yO, t3, x3, y3, t4, x4, y4, v2, False]], Ics = False] ] ; Ics] 
CoordinateChange[tl_, xl_, yl_, t2_, x2_, y2_, tpl_, xpl_, ypl_, tp2_, xp2_, yp2_, 
tp3_, xp3_, yp3_, tp4_, xp4_, yp4_] : = Module [ {atl, all, al2, a21, a22, bl, b2, b3}, 
bl =tl; b2 =xl; b3 = yl; atl = Sqrt[(x2-xl) ^2+ (y2-yl) ^2]; 
all = x2 - xl; al2 = y2 - yl; a21 = yl - y2; a22 = x2 - xl; If [yl == y2, 
f [t_, x_, y_] : = Module [{tt, tx, ty), tt = t; tx = x; ty = y; (tt - bl, tx -b2, ty -b3}] , 
f [t_, x_, y_] : = Module[{tt, tx, ty), tt = t; tx = x; ty = y; 

{atl * (tt - bl) , (tx - b2) * all + al2 * (ty - b3) , a21 * (tx - b2) + a22 * (ty - b3) } ] ] ; 
Join[f[tpl, xpl, ypl], f[tp2, xp2, yp2] , f[tp3, xp3, yp3], f[tp4, xp4, yp4]]] 
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ParamY[x_, t2_, x2_, vl_, t3_, x3_, y3_, v2_] : = 
If[y3j! 0, (((2*x*x2-x2*2 + vl^2*t2''2) / (2*7^2 *t2) -t3) ^2*72^2 - (x-x3)''2- 

vl''2* ((2*x*x2-x2^2 + vl*2*t2*2) / (2 * vl^ 2 *t2) ) -^2 - x''2 - y3''2) / (2*y3), 
Sqrt[(vl* (2*x*x2-x2''2 + vr2*t2''2) / (2 *vl''2 *t2) ) ''2 - x'^2] ] 
ParamYRoot[x_, t2_, x2_, vl_] : = 

(vl* (2* x*x2-x2*2+ 71*2^2-^2) / (2 * vl*2 *t2) ) *2 - x''2 2 
ParamT[x_, t2_, x2_, vl_] := (2 *x*x2 - x2*2 + vl*2 *t2*2) / (2*vl*2*t2) 
ComputeRoots[t2_, x2_, vl_, t3_, x3_, y3_, v2_] : = 

Roots[(2*y3)'^2* ((2*x*x2-x2''2 + vl''2*t2*2)''2*vl*2/ (2*vl*2*t2)*2-x*2) == 
(((2*x*x2-x2*2 + vl^2*t2''2) / (2*vl*2*t2) -t3)^2*v2^2- (x - x3) "2 - y3*2 - 
((2*x*x2-x2''2+vl''2*t2'^2)^2*vl''2/ (2 * vl''2 *t2) -^2 - x''2) ) -^2, x] 
Case3[tl_, xl_, yl_, t2_, x2_, y2_, vl_, t3_, x3_, y3_, t4_, x4_, y4_, v2_] : = 
!£[-■ (vl*v2 = 0), Module [ {Found : False, tls, xls, yls, t2s, x2s, y2s, t3s, x3s, y3s, 
t4s, x4s, y4s}, {tls, xls, yls, t2s, x2s, y2s, t3s, x3s, y3s, t4s, x4s, y4s} = 
CoordinateChange[tl, xl, yl, t2, x2, y2, tl, xl, yl, t2, x2, y2, t3, x3, y3, t4, x4, y4]; 
Iroots : ComputeRoots[t2s, x2s, vl, t3s, x3s, y3s, v2]; 
For[i = 0, (-.Found) A (i< Length [Iroots]), If [lroots[ [i] ] [ [2] ] ifNonreal, 
Found:False, Found:False, Found = ParamYRoot [Iroots [[i]] [[2]], t2s, x2s, vl] A 
(OiParamT[lroots[[i]][[2]], t2s, x2s, vl] i t2s) A 

(CheckHalfSpace[PararaT[lroots[[i]] [[2]], t2s, x2s, vl], Iroots [[i]] [[2]], ParainY[ 
lroots[[i]] [[2]], t2s, x2s, vl, t3s, x3s, y3s, v2], t3s, x3s, y3s, t4s, x4s, y4s, 
v2. True] VIf[y3s " 0, (0 i ParamT[lroots[ [i] ] [ [2] ] , t2s, x2s, vl] S t2s) A 
CheckHalfSpace[ParamT[lroots[[i]] [[2]], t2s, x2s, vl], Iroots [ [i] ][ [2] ] , 
-ParamY[lroots[[i]] [[2]], t2s, x2s, vl, t3s, x3s, y3s, v2], 
t3s, x3s, y3s, t4s, x4s, y4s, v2. True], False])], i++]; 
If [n Found, Module[{}, Iroots = ComputeRoots[t2s, x2s, vl, t4s, x4s, y4s, v2]; 
For[i:0, (-.Found) A (i<Length[lroots]), If [lroots[[i]] [[2]] jiNonreal, 
Found: False, Found: False, Found = ParamYRoot [Iroots [[i]] [[2]], t2s, x2s, vl] A 
(0 < ParamT [Iroots [ [i] ] [ [2] ] , t2s, x2s, vl] < t2s) A 

(CheckHalfSpace[ParamT[lroots[[i]] [[2]], t2s, x2s, vl], Iroots [ [i] ][ [2] ] , 
ParamY[lroots[[i]] [[2]], t2s, x2s, vl, t3s, x3s, y3s, v2], 
t3s, x3s, y3s, t4s, x4s, y4s, v2. False] Vlf [y4s ~ 0, 
CheckHalfSpace[ParamT[lroots[[i]] [[2]], t2s, x2s, vl], Iroots [[i]] [[2]], 
-ParamY[lroots[[i]] [[2]], t2s, x2s, vl, t3s, x3s, y3s, v2], t3s, 
x3s, y3s, t4s, x4s, y4s, v2. False], False])], i++]]. Found = True]; 
If [-.Found, {t3s, x3s, y3s, t4s, x4s, y4s, tls, xls, yls, t2s, x2s, y2s} = 
CoordinateChange[t3, x3, y3, t4, x4, y4, t3, x3, y3, 
t4, x4, y4, tl, xl, yl, t2, x2, y2]. Found: True]; 
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If Found, Module[{}, Iroots = ComputeRoots[t4s, x4s, v2, tls, xls, yls, vl]; 
For[i = 0, (-Found) A (i<Length[lroots]), If[lroots[[i]] [[2]] ^Nonreal, 
Found = False, Found = False, Found = ParamYRoot [Iroots [ [i] ] [ [2] ] , t4s, x4s, v2] A 
(0 iParamT [Iroots [[i]] [[2]], t4s, x4s, v2] i t4s) A 

(CheckHalfSpace[ParamT[lroots[[i]] [[2]], t4s, x4s, v2], Iroots [ [i] ][ [2] ] , 
ParamY[lroots[[i]] [[2]], t4s, x4s, v2, tls, xls, yls, vl], tls, xls, yls, t2s, 
x2s, y2s, vl. True] Vlf [yls == 0, CheckHalfSpace[ParamT[lroots[ [i] ] [ [2] ] , t4s, 
x4s, v2], lroots[[i]] [[2]], -ParamY [Iroots [ [i] ][ [2] ] , t4s, x4s, v2, tls, 
xls, yls, vl], tls, xls, yls, t2s, x2s, y2s, vl. True], False])], i++]; 
If [-■ Found, Module[{}, Iroots = ComputeRoots[t4s, x4s, v2, t2s, x2s, y2s, vl]; 
For[i = 0, (-Found) A (i< Length [Iroots]), If [lroots[ [i] ] [ [2] ] jSNonreal, 
Found = False, Found = False, Found = ParamYRoot [Iroots [ [i] ] [ [2] ] , t4s, x4s, v2] A 
(0iParamT[lroots[[i]][[2]], t4s, x4s, v2] S t4s) A 

(CheckHalfSpace[ParamT[lroots[[i]] [[2]], t4s, x4s, v2], Iroots [[i]] [[2]], 
ParamY[lroots[[i]] [[2]], t4s, x4s, v2, t2s, x2s, y2s, vl], tls, xls, yls, t2s, 
x2s, y2s, vl. False] VIf[y2s==0, CheckHalfSpace [ParamT [Iroots [ [i] ] [[2]], 
t4s, x4s, v2], Iroots [[i]] [[2]], -ParamY[lroots[ [i] ] [ [2] ] , t4s, x4s, v2, 
t2s, x2s, y2s, vl], tls, xls, yls, t2s, x2s, y2s, vl. False], False])], 
i++] ] , Found = True] ] , Found = True] ; Found] ' 
alibi [tl_, xl_, yl_, t2_, x2_, y2_, vl_, t3_, x3_, y3_, t4_, x4_, y4_, v2_] : = 
If [isBeadNonEmpty[tl, xl, yl, t2, x2, y2, vl] A isBeadNonEmpty[t3, x3, y3, t4, x4, y4, v2], 
If[Casel[tl, xl, yl, t2, x2, y2, vl, t3, x3, y3, t4, x4, y4, v2]. True, 
If[Case2[tl, xl, yl, t2, x2, y2, vl, t3, x3, y3, t4, x4, y4, v2]. True, 
Case3[tl, xl, yl, t2, x2, y2, vl, t3, x3, y3, t4, x4, y4, v2]]]. False] 
alibiQE[tl_, xl_, yl_, t2_, x2_, y2_, vl_, t3_, x3_, y3_, t4_, x4_, y4_, v2_] : = 
Resolve [Exists [{t, x, y}, 

((x-xl)^2+ (y-yl)*2ivl*2*(t-tl)*2A (x-x2)*2+ (y - y2) *2 5 vl*2 * (t - 12) *2 A 
(tl it i t2)) A ((x-x3)^2+ (y-y3) "2 i v2^2 * (t - 13) ^2 A 
(x-x4)*2+ (y-y4) ^2 i v2*2* (t-t4) '^2A (t3it2t4))]] 
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